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PREFACE 


Undaratanding  tha  various  regions  of  •  glow  diacharga  haa 
lad  to  laprovaaanta  in  a  wida  variaty  of  military  ayatama.  auch 
aa  high  powar  axciaar  laaara,  thyratrona  for  high  powar  awitchaa, 
and  plaama  raactora  for  atching  or  dapoaition  of  thin  filma  in 
t.ha  aanufactura  of  aaalconductor  ehipa.  Many  of  thaaa  diachargaa 
of  currant  military  intaraat  now  contain  raactiva  alactronagativa 
gaaaa.  It  ia  hopad  th 1 a  manuacrlpt  haa  providad  tha  raadar  a 
battar  undaratanding  of  tha  influanca  alactronagativa  gaaaa  hava 
on  tha  cathoda  fall  ragion  of  a  glow  diaeharga.  Thia  manuacnpt 
aumaarizaa  tha  raaulta  of  a  mora  coaprahanaiva  AFIT  doctoral 
diaaartation . 

1  grataful  to  Or.  Alan  Garacaddan.  my  diaaartation 
advisor  at  tha  Aaro  Propul axon  Laboratory,  and  Haior  Douglaa 
Daablar,  Air  Command  and  Staff  Collaga  faculty  mambar ,  for 
rssdlng  thia  manuacript  and  providing  thair  critical  commanta. 
Spseisl  thanka  goaa  to  tha  Aaro  Propulaion  Laboratory  Powar 
Diviaion  who  aponaorad  thia  atudy. 

Sub? act  to  clasrsnca.  thia  aanuacript  will  ba  aubaittad  to 
tha  JgyEnii_gf_^ppiigd_Phy§lg§  for  conaidarat ion .  Tha  format  of 
tha  manuacript.  including  rafarancaa,  ara  conaiatant  with  tha 
raquiraaanta  for  publication  in  that  journal . 
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EXECUTIVE  SUMMARY 


"“insights  into  tomorrow ” 


Part  of  our  College  mission  is  distribution  of  the 
students’  problem  solving  products  to  DoD 
sponsors  and  other  interested  agencies  to 
enhance  insight  into  contemporary,  defense 
related  issues.  While  the  College  has  accepted  this 
product  as  meeting  academic  requirements  for 
graduation,  the  views  and  opinions  expressed  or 
implied  are  solely  those  of  the  author  and  should 
not  be  construed  as  carrying  official  sanction. 


REPORT  NUMBER  67-0730 

AUTHOR(S)  MAJOR  GARY  L.  DUKE.  PhD.  USAF 

TITLE  INVESTIGATION  OF  SHEATH  PHENOMENA 

IN  ELECTRONEGATIVE  GLOW  DISCHARGES 

I.  Pu rgoae :  To  educate  the  scientific  community  about  now  the 
addition  of  a  qaa  wtucft  easily  torma  negative  ions  'electronega¬ 
tive  gas)  influences  the  cathode  region  and  thus  the  voitaoe  ciron 
across  a  glow  discharge.  A  glow  discharge  consists  ot  a  set  or 
electrodes  surrounded  by  a  gas  or  gas  mixture.  A  fluorescent 
light  is  a  simple  glow  discharge  operated  with  an  alternatino 
current . 

II.  Problem:  The  width  of  the  iight  and  dark  reaions  ciose  to 
the  cathode  have  been  experimentally  observed  (Rei  &>  to  contract 
longitudinally  toward  the  cathode  when  an  electronegative  aas  was 
added  to  a  rare  (inert)  gas  discnarae.  A  smaller  voitaoe  drop 
across  this  region  is  anticipated  as  a  result  of  this 
contraction.  Theoretical  models  constructed  to  describe  electric 
discharge  lasers,  high  power  switches,  and  devices  used  to 
manufacture  semiconductors  have  typically  treated  tne  voitaoe 
drop  across  the  cathode  fall  region  as  a  constant,  as  if  the 
electronegative  gas  did  not  affect  it.  This  paper  attempts  to 
resolve  this  inconsistency  by  explaminq  why  the  width  of  tne 
cathode  fall  contracts  when  an  electronegative  aas  is  aooed  to 
the  discharge. 

III.  Data :  A  theoretical  model  of  the  cathode  fall  region  was 
constructed.  The  electron  kinetics  were  analyzed  by  soivma  the 
one  dimensional,  time  independent,  collisional  Boltzmann  equation 


CONTINUED. 


which  describes  statistically  a  swarm  of  alectrons  as  it 
traverses  tha  cathode  fall  region.  Tha  only  production  and  loss 
machanisms  for  alactrons  was  assumad  to  ba  ionization  and 
attachmant.  No  loss  machanism  was  assumad  to  exist  for  negative 
ions.  This  assumption  holds  for  low  pressure  discharges  where 
ion-ion  recombination  can  ba  neglected.  Tha  positive  ion 
kinetics  are  based  on  Poisson" s  aquation  which  relates  the  net 
charge  density  to  a  fluctuation  in  tha  electric  field  and  an 
approximation  for  tha  mobility  of  tha  ions  in  high  electric 
fields.  Tha  negative  ion  density  is  derived  from  the  electron 
kinetics,  since  tha  negative  ions  are  formed  by  slow  electron 
collisions  with  tha  electronegative  gas  molecules.  The  negative 
ion  mobility  is  assumad  to  ba  the  same  as  the  positive  ion 
mobility.  This  is  a  good  first  approximation,  but  could  be 
refined  in  a  later  analysis. 

Gas  mixtures  of  helium,  argon,  and  xenon  with  lx  and  5x 
hydrochloric  acid  were  analyzed.  Tha  electric  field,  the 
Townsend  ionization  and  attachment  rates,  and  charged  particle 
densities  were  calculated  as  a  function  of  distance  across  the 
cathode  fall  region  and  compared.  A  contraction  was  predicted  in 
the  helium  mixtures,  none  was  predicted  in  the  argon  mixtures, 
and  a  slight  expansion  in  the  cathode  fall  was  predicted  in  xenon 
mixtures. 

IV.  I  Tha  results  indicated  ionization  and  not 

attachment  of  a  gaa  auch  as  hydrochloric  acid  m  helium  mixtures 
led  to  the  contraction  and  a  reduced  voltage  drop  across  the 
cathode  fall.  Since  the  ionization  threshold  of  the  added 
electronegative  gas  was  less  than  the  rare  gas,  less  voltage  was 
required  to  accelerate  electrons  to  achieve  the  same  net 
ionization  in  this  region.  It  was  also  observed  that  the  gas 
with  the  lower  ionization  threshold  was  spatially  ionized  first 
leading  to  a  separation  of  the  gas  components  in  the  cathode 
region.  The  density  of  positive  ions  resulting  from  electron 
ionization  at  the  edge  of  the  cathode  fall  was  always  found  to 
determine  the  slope  of  the  electric  field  in  all  mixtures. 

Conversely,  results  indicated  electron  attachment  in  the 
cathode  fall  region  increased  the  cathode  fall  voltage  and 
length.  Although  the  attachment  coefficient  was  larger  than  the 
ionization  coefficient  immediately  in  front  of  the  cathode,  the 


viii 


CONTINUED 


density  of  negative  Iona  did  not  contribute  aignificantly  to  the 
apace  charge  until  near  the  end  of  the  cathode  fall.  The  voltage 
drop  and  length  of  the  cathode  had  to  increaae  in  order  to  main¬ 
tain  the  aame  net  ionization  aa  occurred  in  a  pure  rare  gaa.  The 
claaaical  theory  that  the  formation  of  negative  iona  in  the 
cathode  fall  cauaea  it  to  contract  due  to  their  slower  mobility 
waa  shown  to  be  falae. 

V.  Recommendatlona :  Scientists  and  engmeera  modeling  electric 
diacharge  laaera,  high  power  awltchea,  and  plaama  deposition  and 
etching  devicea  ahould  not  aaaume  the  voltage  drop  across  the 
cathode  fall  region  in  discharges  containing  mixtures  is  equiva¬ 
lent  to  the  voltage  drop  acroaa  the  cathode  fall  region  of  a 
diacharge  containing  the  primary  conatituent  (usually  the  rare 
gaa) .  The  voltage  drop  acroaa  the  cathode  fall  can  be 
aignificantly  affected'by  trace  species  if  their  threshold  for 
ionization  is  leas  than  the  primary  conatituent.  Scientists 
modeling  plaama  deposition  and  etching  devicea  should  also  be 
aware  that  the  gaaea  are  sequentially  ionized  across  the  cathode 
fall  region  in  order  of  their  ionization  thresholds.  This 
phenomena  could  impact  the  uniformity  and  reproducibility  of 
their  films. 


IX 


I .  INTRODUCTION 


Previous  papers  have  analyzed  -the  cathode  fall  in  electro¬ 
positive  glow  discharges.!’  2*  3»  4»  5*  6»  7  This  paper  analyzes 
the  cathode  fall  in  an  electronegative  glow  discharge.  An 
understanding  of  the  formation  of  positive  and  negative  ions  in 
electronegative  glow  discharges  is  very  important  in  modeling 
electric  discharge  lasers,  high  power  switches,  and  the  plasma  in 
plasma  etchers  and  reactors  for  depositing  thin  films.  In  1938, 
Emeleus  and  Sayers®  reported  the  width  of  the  luminous  and  dark 
regions  in  the  cathode  fall  and  negative  glow  region  of  a  helium 
discharge  contracted  significantly  towards  the  cathode  with  the 
addition  of  just  a  trace  of  chlorine,  an  electronegative  gas. 

This  contraction  of  the  cathode  fall  implies  a  decrease  m 
cathode  fall  voltage.  This  change  in  voltage  is  important 
because  it  determines  the  energy  of  the  electrons  and  ions  which 
traverse  the  cathode  fall  region  and  enter  the  negative  glow. 

This  paper  explains  the  contracted  cathode  fall  region  observed 
in  electronegative  glow  discharges. 

Two  theoretical  approaches  exist  for  calculating  the 
nonequilibrium  electron  energy  distribution  function  in  the 
cathode  fall  region  of  a  glow  discharge.  The  Monte  Carlo  method 


is  based  upon  the  statistics  of  modeling  a  large  number 


(typically  10^  or  greater)  of  electron  collisions  with  gas  mole¬ 
cules  in  the  presence  of  an  electric  field. 3*  7 ’  9  The  other 
method  is  based  upon  numerical  1-y  solving  the  time  independent 
Boltzmann  equation  for  the  electron  distribution  function.  In 
uniform  electric  fields,  agreement  between  these  methods  has  been 
excellent.^0'  !!'  12  Allis  and  co-workers! '  2  and  Long^  assumed 
a  linearly  decreasing  electric  field  through  the  cathode  fall  and 
used  the  Boltzmann  method  to  model  the  change  in  the  electron 
distribution  function  as  a  function  of  distance.  Tran  Ngoc,  et 
al.?  and  Segur^  similarly  assumed  a  linearly  decreasing  electric 
field,  but  used  the  Monte  Carlo  model.  These  calculations  are 
not  considered  self-consistent.  The  perturbation  of  the  electric 
field  from  electron  nonequilibrium  kinetics,  which  their 
calculations  showed  had  occurred  in  the  cathode  fall  region,  was 
not  examined.  Since  then  Segur  et  al.6  have  coupled  a  Boltzmann 
solution  for  electrons  with  an  analytical  approximation  to  the 
Boltzmann  equation  for  positive  ions  in  the  cathode  fall  region 
in  a  self-consistent  manner  for  various  helium-mercury  mixtures. 
Boeuf  et  al.4  investigated  the  cathode  fall  in  a  helium  discharge 
using  a  self-consistent  Monte  Carlo  technique,  which  sequentially 
iterated  Monte  Carlo  calculations  for  positive  ions  followed  by  a 
calculation  of  the  electric  field,  and  then  another  Monte  Carlo 
calculation  for  electrons.  The  present  paper  is  based  on  the- 
Boltzmann  technique  derived  by  Allis  et  al.2  and  Long^,  and 
modified  to  include  the  formation  of  negative  ions  as  well  as  a 
self-consistent  calculation  of  the  electric  field. 
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II 


MODELING 


a .  Thaor y 

Tha  cathode  fall  region  ia  assumed  to  be  a  weakly  ionized, 
collision  dominated,  nonneutral  plaaaa  consisting  of  four  compon- 
enta:  electrons,  positive  ions,  negative  ions,  and  neutraia. 

The  densities  of  negative  and  positive  ions  are  derived  from  the 
negative  ion  current  density  and  current  conservation  respec¬ 
tively.  The  neutral  density  is  assumed  to  remain  constant. 

The  electron  density  and  the  electron  energy  distribution 
are  determined  by  numerically  solving  the  time  independent,  one 
dimensional,  colliaional  Boltzmann  equation.  This  representation 
models  the  distribution  f  of  a  swarm  of  electrons  as  it  traverses 
the  cathode  fall  of  a  glow  discharge.  The  one  dimensional 
Boltzmann  equation 


is  transformed  into  energy  apace  as 


<1> 


1L  +  11  =  l-  if  ) 
\P  )S  Jt  I,, 


<2) 


where  0*ejE<x)dx  is  the  electron  potential  energy  and  £  =»vx2/2  is 
the  electron  kinetic  energy  induced  by  the  electric  field  E<x>. 
Let  the  total  electron  kinetic  energy  be  represented  by  £=£•►? 


wh«r«  t  la  th«  random  or  tharmal  anargy  cont.ribut.ion.  Tha 
diatribution  function  f  can  than  ba  rapraaantad  aa  tha  aum  of  two 
functions,  f-<£,£,0),  which  ia  tha  diatribution  of  alactrona  with 
vx  < O  moving  againat  tha  fiald  and  f+<£,?,0>,  which  rapraaanta 
thoaa  alactrona  with  vx > O  moving  with  tha  fiald.  Tha  right  hand 
tarm  rapraaanta  tha  changa  in  f  dua  to  tha  colliaion  procaaaaa  in 
tha  gaa.  Thia  aquation  aubaaquantly  ia  tranaformad  into  a  linaar 
intagral  oparator  aquation.13 


whara  K  ia  an  intagral  oparator  rapraaanting  a  aum  ovar  all 
colliaional  procaaaaa,  Q<1>  ia  tha  total  alactron  croaa  aaction, 
and  N  ia  tha  gaa  danaity.  Darivation  of  tha  individual  oparatora 
for  aach  alactron  colliaional  procaaa  can  ba  found  in  Long' a 
raport . 5 

Tha  colliaional  tarm  ia  furthar  aimplifiad  by  aaauming  Kf  to 
ba  naarly  iaotropic  in  valocity  apaca  and  axpanding  Kf  in 
Lagandra  polynomiala.  Thia  axpanaion  diffara  from  tha  uaual  two 
tarm  axpanaion1  of  Boltzmann' a  aquation  in  apharical  harmonica  in 
that  it  ia  an  axpanaion  of  tha  product  of  tha  intagral  oparator 
timaa  f ,  not  an  axpanaion  of  f.  Tha  raaulting  aquation  ia  aolvad 
numarically  uaingan  itarativa  tachniqua  davalopad  by  Long5.  It 
ia  aaaumad  initially  that  alactrona  are  only  scattarad-out  of 
aach  alamant  of  valocity  apaca  <Kf *0> .  Tha  aquation  ia  than 
intagratad  for  both  forward  and  backward  moving  alactrona. 
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Scattered -out  electrons  than  becoae  scattered-in  alactrons  for 
othar  alasants  of  valocity  apaca.  Thia  givaa  a  firat  approxima¬ 
tion-  to  tha  alactron  distribution  function  fros  which  a  naw  oper¬ 
ator  Kf  ia  calculated.  Thia  procadura  waa  rapaatad  until  conver¬ 
gence  waa  reached  (taken  aa  laaa  than  . 2*  change  in  tha  diatnbu- 
tion  function).  For  a  200  aV  cathode  fall,  tha  energy  raaolution 
of  tha  calculation  corraapondad  approximately  to  1  eV . 

From  tha  alactron  energy  diatr ibut ion ,  tha  alactron  number 
denaity ,  alactron  currant  denaity,  and  Townaand  ionization 
coefficient  ware  calculated.  Tha  alactron  energy  and  number 
denaity  diatr ibutiona  ware  normalized  to  tha  diatnbution  at  tha 
cathode.  In  order  to  obtain  a  completely  aalf -conaiatant 
aolution  in  which  Poiaaon'a  aquation  la  obeyed  in  addition  to 
boundary  conditiona  at  tha  cathode  and  anode,  successive  calcula¬ 
tions  for  tha  alactron  diatnbution  function,  electron  number 
denaity,  and  raaulting  electric  field  ware  accompl iahed . 

The  electric  field  ia  obtained  by  putting  Poisson's  equation 
in  difference  form  ao  tha  value  of  the  field  at  each  point  is 
baaed  on  the  value  at  the  previous  increment  with  the  field  at 
the  cathode  given  aa  a  boundary  condition. 

Ej  -  Ej.(  ~  ”  *V/  “  -  fi-J  “  n-iv 

where  n*j,  n«j ,  and  n-j  are  the  positive  ion,  electron,  and 
negative  ion  densities  at  position  i,  and  xi  is  the  spatial  point 
corresponding  to  increment  j.  The  positive  ion  density  is 


written  in  taraa  of  tho  currant  conservation  (J  ■  total  currant 
danaity)  using  Ward's  sodifisd  formula14  for  tho  ion  drift 
velocity  for  high  field  regions. 


where  3*3  and  3-3  are  the  electron  and  negative  ion  current 
densities  at  location  3.  Similarly  the  negative  ion  density  is 
calculated  by  dividing  the  negative  ion  current  density  oy  its 
velocity  ( k -  is  assumed  to  be  equal  to  k*>. 


J- 


Since  the  negative  ions  are  formed  as  a  result  of  electron 
attachment,  the  negative  ion  current  density  is  obtained  oy 
integrating  over  all  attaching  collisions. 

« 

1  *  J/3,  J.t  A 

where  1a  the  attachment  rate.  Equation  4  then  becomes  a  fifth 
order  polynomial  whose  roots  were  found  numerically  using 
Steinman'a  method. ^5  The  lowest  real  value  of  E  (some  are 
imaginary)  at  the  cathode  boundary  was  picked  as  the  solution, 
since  it  led  to  the  lowest  possible  cathode  tall  voitage.  The 
electric  field  at  each  location  is  therefore  calculated  from  ne. 
3«,  3-  and  the  value  of  the  electric  field  at  the  previous 


•patlal  liicrMiiit .  This  tschnlqus  eenvsrysd  quickly,  neraslly 
within  4-4  itsrstlocia  for  thrss  olgnlficent  flgura  convargsncs . 

b.  Coo  and  D lecher go  Psrssotsrs 

Tteo  gas  and  diaehsrgs  psrasstsrs  uood  in  tho  colculotiono 
oro  suassrlasd  in  Toblo  1 .  For  •  roquirod  sloctron  energy 
rooolution  of  1  eV,  tho  oloetroda  dlotoneo  (D>  woo  linitod  by  tho 
ooo or y  etorege  copobility  of  o  Cybor  74  coaputar.  For  100*  Ho,  D 
woa  2.2  ea  ond  for  tho  roooininq  Ho/HCl  aixturoo,  it  woo  2.1  ca 
duo  to  tho  eontroction  thot  occurred .  t  io  tho  eecondery 
electron  ooisaion  coefficient  of  tho  cothodo  for  ion  boeberdeent. 
Tho  dietribution  of  oloetrono  looving  tho  cothodo  woo  tho  eeee  aa 
thot  used  by  Long9  (14*  between  1-2  eV,  25*  between  2-3  eV,  24* 
between  3-4  oV,  25*  between  4-5  oV,  ond  4*  between  5-4  oV> .  Tho 
nobility,  k* ,  of  tho  poeitlve  iono  woo  repreeented  by  those  of 
tho  roro  gone*. 

Toblo  I.  Coo  ond  dlecherge  poroootoro  uood  in  tho  Boltxaonn 
colculotiono. 


Mixture* 

D 

(ca) 

J 

(aoap/ca2) 

r 

k 

ca-torr 

volt-ooc 

Ho/HCl 

2. 1-2. 2 

14 

.2 

4. IE-4 

Ar/HCl 

1.15 

10 

.0417 

4 . 25E ♦ 3 

Xo/HCl 

.47 

10 

.004 

4 . OE*  3 

c.  Electron  Zopoct  Crooa-Soctiona. 

Tho  cross- sect ions  uood  for  Ho,  Ar ,  Xo,  ond  HC1  in  tho 
boltzoonn  colculotiono  oro  shown  in  Figo.  1-5.  Threaholdo  and 
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Figure  3  ELECTRON  IMPACT  CROSS  SECTIONS  IN  Xe. 
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Figure  3  LOW  ENERGY  ELECTRON  CROSS  SECTIONS  IN  HCI. 


rafarancaa  for  aach  anargy  rang*  ara  liatad  in  Tabla  II.  Tha 
aoat  important  croaa  aactiona  in  tha  calculationa  ara  tha 
ionization  croaa  aactiona  of  tha  rara  gaaaa.  Thia  croaa  aaction 
datarainaa  how  quickly  tha  alactrona  multiply*  and  thua  diractly 
af facta  n«  and  dE/dx  in  Poiaaon'a  aquation.  For  tha  rara  gaaaa* 
tha  ionization  croaa  aaction  waa  obtainad  from  maaauramanta  by 
Rapp  at  al.16  ovar  tha  complata  anargy  ranga.  Tha  ionization 
croaa  aaction  for  HC1  had  to  ba  axtandad  abova  lOO  aV.  Sinca  HC1 
waa  uaad  in  amall  concantrationa  and  tha  maan  alactron  anargy  in 
tha  HC1  mixturaa  waa  alwaya  laaa  than  87  aV*  any  arror  in  tha 
approximation  ahould  ba  nagligibla.  Elactronic  croaa  aactiona 
for  Ha  wara  axtandad  uaing  tha  formula*  liatad  undar  tha 
rafaranca  column  in  Tabla  II. 


III.  DISCUSSION 


Tha  thaory  daacribad  abova  waa  uaad  to  modal  tha  cathoda 
fall  in  Ha/HCl *  Ar/HCl,  and  Xa/HCl  mixturaa.  All  data  praaantad 
waa  obtainad  for  a  gaa  praaaura  of  1  torr  at  a  tamparatur*  of  300 
K.  Gaa  and  diacharg*  paramatara  pravioualy  daacribad  in  Tabla  I 
war*  hald  conatant.  Tha  only  production  and  loaa  procaaaaa  for 
alactrona  waa  aaauaad  to  b*  ionization  and  attachmant.  At  low 
gaa  praaauraa  ion-ion  racombination  can  b*  naglactad.  No  loaa 
procaaa  for  nagativa  iona  waa  aaauaad  to  axiat. 


Table  II.  Thresholds  for  electron  impact,  cross  sections. 


1  Procsss 

Range  <  e V  > 

References 

Hs 

Hossntus 

0.-6. 

17 

Transfer 

6.-500. 

18 

Electronic 

23s  ♦  2lS 

19.81-23.54 

17 

24.-100. 

Constant 

100.-500. 

«■  ■  .624/6* 

23p 

20.95-199.6 

17 

200.-500. 

C*  »  .02755/6*- 

2lp 

21.2-196.1 

17 

200.-500. 

<r  *  1.179/6* 

3lp 

23.07-207.8 

17 

210.-500. 

er  »  .3605/6* 

Ionization 

24.59-500. 

16 

Ar 

Moaantus 

0.-30. 

17 

Transfer 

30.-500. 

5 

Electronic 

2**3/2<3PS> 

11.6-500. 

19 

2pi/2<2p5) 

11.8-500. 

19 

Remaining 

13.2-500. 

19 

Ionization 

15.76-500. 

16 

Xe 

Momentum 

0.-6.5 

17 

Transfer 

6.5-20. 

19 

20.-500. 

18 

Electronic 

8.32-500. 

19 

Ionization 

12.13-500. 

16 

HC1 

Momentum 

0.-100. 

20 

Tr ansf er 

100.-500. 

Log  Extrapolation 

Vibration 
v  »  1 

K 

i 

si 

. 

20 

v  »  2 

.7-7.2 

20 

Attachment 

Cl- 

.67-2.7 

20 

H* 

5.6-12. 

20 

Excitation 

A 

5.5-11. 

20 

B  ♦  C 

9.3-100. 

20 

100.-500. 

19 

Ionization 

12.74-100. 

20 

100. -500. 

19 

14 


a.  Helium  aixturaa 


A  15*  dacraaaa  in  the  a lac trie  fiald  ia  pradictad  aa  tha 
amount.  o£  HC1  ia  incraaaad  froa  0  to  5*  in  Fig.  6.  Thia 
corraaponda  to  a  32*  dacraaaa  in  cathoda  fall  voltaga  and  a  21* 
dacraaaa  in  cathoda  fall  langth.  Thia  dacraaaa  in  voltaga  and 
diatanca  aa  a  function  of  HC1  ia  dua  to  tha  incraaaa  in  tha 
Townaand  ionization  coefficient,  aa  diaplayad  in  Fig.  7.  Thia 
incraaaa  in  ionization  in  tha  paak  ragion  occura  bacauaa  HC1  ia 
much  aaaiar  to  ioniza  than  Ha  (tha  ionization  thraahold  for  HC1 
ia  12.74  aV ,  varaua  24.59  aV  for  Ha) .  Also,  tha  apatial 
thraahold  for  tha  Townaand  ionization  coafficiant  ia  halvad  whan 
HC1  ia  added,  raflacting  tha  changa  in  thraahold  for  ionization. 
Thua  vary  cloaa  to  tha  cathoda,  aora  HCL  ia  baing  ionizad  than 
Ha.  Thia  ia  conaiatant  with  tha  raaulta  of  Segur,  at  al.6  with 
Hg  in  Ha. 

Tha  alopa  of  tha  alactric  fiald  naar  tha  cathoda  ia 
daacribad  by  Poiason's  aquation.  Aa  will  ba  aaan  in  latar 
figuraa,  tha  poaitiva  ion  danaity  naar  tha  cathoda  ia  auch  larger 
than  othar  charged  particle  denaitiaa  by  aaveral  ordera  of 
magnitude.  The  linear  alopa  of  the  electric  field  thua  dependa 
only  on  the  near  conatant  value  of  the  poaitive  ion  danaity.  The 
curvaa  in  Fig.  6  are  thua  almoat  parallel.  The  electric  field 
curvea  are  ahiftad  down  becauae  leaa  voltage  ia  required  to 
produce  the  aama  net  ionization. 

Townaand  ionization  and  attachment  coefficienta  are  compared 
in  Fig.  8  aa  a  function  of  diatanca  for  99/1  and  95/5  concentra- 
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ELECTRIC  FIELD  IN  He/HCl  MIXTURES. 


DISTANCE  (cm) 

I  ATTACHMENT  COEFFICIENTS  IN  He/HCI  MIXTURES. 


tiona  of  He/HCl 


At  ont  would  txptet,  tht  ttgnltudt  of  tht 


attachment  rata  la  auch  aaallar  than  tha  ionization  rata 
throughout  tha  cathode  fall  region,  axcapt  iaaadiataly  in  front 
of  tha  cathoda.  In  thia  ragion,  which  corraaponda  to  tha  Faraday 
dark  apace,  alactrona  hava  juat  laft  tha  cathoda  with  a  faw  eV  of 
anargy  and  ara  aora  auacaptibla  to  attachaant  than  aftar  thay 
hava  baan  accalaratad  to  highar  anargiaa.  Tha  aacond  paak  aaan 
in  tha  attachaant  curva  at  1  ca  ia  probably  dua  to  tha  alactron 
diatribution  function  firat  raaching  tha  ionization  anargy  and 
producing  a  naw  crop  of  alow  aacondary  alactrona.  Thia  ia 
conaiatant  with  tha  paak  in  tha  ionization  coafficiant  which 
occura  naar  tha  aaaa  location. 

Tha  alactron,  poaitiva  ion,  and  nagativa  ion  nuabar 
danaitiaa  ara  dlaplayad  aa  a  function  of  diatanca  through  tha 
cathoda  fall  ragion  in  Fig.  9.  Tha  alactron  danaity  initially 
dacraaaaa  bacauaa  tha  majority  of  alactrona  ara  rapidly 
accalaratad  away  from  tha  cathoda.  Juat  bafora  thay  hava  gainad 
aufficiant  anargy  for  ionization,  thair  danaity  raachaa  a 
minimum.  Tha  alactron  danaity  bagina  to  incraaaa  aftar  tha 
ionization  potantial  of  tha  gaa  ia  attainad.  Similarly,  tha 
poaitiva  ion  danaity  incraaaaa  alightly  laaving  the  cathoda 
bacauaa  vary  littla  ionization  occura  cloaa  to  tha  cathoda  aa  tha 
poaitiva  iona  ara  baing  accalaratad  towarda  it.  Although 
nagativa  iona  ara  praaant  throughout  tha  cathode  fall  ragion, 
thair  danaity  ia  ganarally  much  too  amall  to  affect  the  net  apace 
charge  to  any  aignlf leant  degree.  Even  in  tha  95/5  mixture  tha 
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Figure  *  C  HARGED  PARTICLE  DENSITIES  IN  A  *3/5  He/HCI  MIXTURE. 


Mfatlv*  ion  tensity  hm  atlll  nor*  thon  on  or  dor  of  negni todo 
looo  thon  tho  oloetron  dowel ty  ond  is  nor*  thon  four  or dors  of 
nnynltute  onollor  thon  tho  positive  ion  density.  Tho  negnitude 
of  tho  oloetron  donoity  opprooehoo  tho  pooitivo  ion  denolty  only 
noor  tho  eothodo  f oil -nogotivo  glow  boundory .  Cguilibriun  woo 
novor  rooehod  in  thooo  onoaploo  boeouoo  o  ought ly  blnodol 
oloetron  diotribution  otill  existed  ot  thin  boundory. 

b.  Argon  nixtures 

In  controot  to  Ho  nixtures,  thoro  lo  o  nogligiblo  offoet  on 
tho  oloctrie  fiold  in  tho  eothodo  foil  roglon  when  up  to  So  HC1 
io  oddod  to  Ar  oo  shown  in  Fig.  10.  Thlo  oloctrie  fiold  diotri¬ 
bution  through  tho  eothodo  foil  corroopondo  to  obout  o  144  v  po- 
tontiol  drop  oerooo  tho  eothodo  foil.  With  1*  HCl,  tho  differ - 
once  in  totol  ionlsotion  io  too  onoll  to  hove  ony  oboorvobio 
offoet  on  tho  oloetric  field.  Ionlsotion  throoholdo  for  Ar  <15.7 
oV)  ond  HCl  (12.7  oV)  ore  ouch  elooor  thon  for  Ho  ond  HCl. 

Thio  difference  io  oloo  too  onoll  to  bo  oboorvod  in  Fig.  11 
in  tho  Townoend  ionlsotion  coefficient.  With  5*  HCl,  thoro  is 
oufflciont  HCl  to  incroooo  tho  ionlsotion  roto  obout  l.€»* 
throughout  ooot  of  tho  eothodo  foil  region.  Townoend  ionlsotion 
ond  ettechnent  coefflclente  for  99/1  end  95/5  concontrot iono  in 
Ar/HCl  ore  oloo  coopered  in  thio  figure.  An  in  Ho/HCl  nixtures. 
the  ettechnent  roto  io  lerger  thon  tho  ionlsotion  roto  only  for  a 
few  electron  collleione  clooo  to  tho  eothodo  whore  tho  electron* 
hove  ineufflclent  energy  for  ionlsotion. 
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Tha  mm  nonequilibrium  pkanoMna  that,  occurred  in  nuabar 
danaitiaa  in  Ha  aixturaa  ia  obaerved  in  A r  mlxturea  in  Fig.  12. 
Tha  poaitiva  ion  danaity  ineraaaaa  in  tha  diraction  moving  away 
from  tha  ea thode,  whila  tha  alactron  numbar  danaity  initially 
daeraaaaa .  Tha  nuabar  of  nagativa  Iona  ia  mora  than  two  ordara 
of  magnituda  amallar  than  tha  nuabar  of  poaitiva  iona,  and 
actually  thalr  danaity  doaa  not  bag in  to  grow  axponantially  until 
tha  cathoda  fall  bagina  to  aarga  into  tha  nagativa  glow.  Tha 
poaitiva  ion  danaity  ia  atill  aavarml  ordara  of  magnituda  largar 
than  tha  alactron  and  nagativa  ion  danaitiaa  and  thua  datarminam 
tha  alopa  of  tha  alaetric  fiald. 

c.  Xanon  aixturaa 

In  contract  to  Ha  and  A r  aixturaa,  thara  ia  a  alight 
incraaaa  in  tha  alaetric  fiald  in  tha  cathoda  fall  ragion  whan  up 
to  Si  HC1  ia  addad  to  Xa  aa  ahown  in  Fig.  13.  Slightly  highar 
fialda  ara  raquirad  to  obtain  tha  aama  nat  ionization  naar  tha 
cathoda  dua  to  tha  formation  of  nagativa  iona  vary  cloaa  to  tha 
cathoda  aa  wall  aa  tha  additional  anargy  loaaaa  to  HC1  vibra¬ 
tional  procaaMa.  Xa/HCl  aixturaa  ara  uniqua  in  that  both 
conatituanta  hava  almoat  tha  aama  ionization  potential  (12.74  mV 
for  HC1  and  12.13  aV  for  Xa) .  A  loaa  of  alactrona  cloaa  to  tha 
cathoda  ia  aimilar  to  tha  effect  of  reducing  tha  aecondary 
amlaaion  coefficient  of  tha  cathoda.  A  dacraaaa  in  if  aimilar ly 
raaulta  in  a  highar  cathoda  fall  voltage.  However ,  tha  alopa  of 
tha  electric  fiald  ataya  tha  aama  bacauaa  poaitiva  iona  atill 
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dominate  the  other  number  denaitiee.  Thua  the  trend  of  the 


electric  field  ahifting  vertically  with  a  corresponding  increaaa 
in  voltage  would  be  expected  to  continue  for  larger  concentra- 
tiona  of  HC1  in  Xe.  In  thia  atudy  the  poaitiva  ion  velocity  of 
only  the  rare  gaa  waa  included.  Addition  of  HC1  with  ita  higher 
mobility  would  increaaa  the  average  poaitiva  ion  drift  velocity. 
Thia  would  reault  in  a  shallower  slope  of  the  electric  field  and 
a  further  expansion  in  cathode  fall  length. 

Townsend  ionization  and  attachment  coefficients  are  compared 
in  Fig.  14  for  99/1  and  9S/5  concentrations  in  Xe/HCl .  Again  the 
attachment  coefficient  is  much  smaller  than  the  ionization  rate 
throughout  the  cathode  fall  region  except  immediately  in  front  of 
the  cathode.  Note  that  the  region  where  the  attachment  coeffi¬ 
cient  is  greater  than  the  ionization  rate  is  smaller  in  Xe 
mixtures  than  in  either  He  or  Ar  mixtures  (Figa.  S  and  11  respec¬ 
tively)  .  The  onset  of  ionization  however  still  occurs  cloae  to 
the  second  peak  of  tha  attachment  rate. 

The  nonequilibrium  phenomena  occurring  in  the  previous  rare 
gaa  mixtures  la  also  observable  in  tha  number  densities  in  Xe 
mixtures.  Results  are  illustrated  in  Fig.  15.  The  trends  are 
similar  to  thoae  seen  in  Figa.  9  and  12. 


IV.  SUMMARY  AND  CONCLUSIONS 

These  results  indicate  ionization  and  not  attachment  of  a 
gaa  such  aa  HC1  in  He  mixtures  leads  to  a  contraction  and  a 
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raducad  voltage  drop  across  ths  cathode  fall.  Since  the 
ionization  threshold  of  the  added  electronegative  gas  was  less 
than  the  rare  gas,  less  voltage  was  required  to  accelerate 
electrons  to  achieve  the  same  net  ionization  in  the  this  region. 
This  explains  the  contraction  of  the  cathode  fall  region  observed 
by  Eaeleus  and  Sayers^  in  electronegative  gas  Mixtures.  It  was 
also  observed  that  the  gas  with  the  lower  ionization  threshold 
was  ionized  first.  This  leads  to  a  separation  of  the  gas 
coeponents  in  the  cathode  fall  region.  Discharge  models  which 
include  the  cathode  fall  voltage  in  their  circuit  analysis  should 
include  the  different  ionization  rates  and  distinct  mobilities 
for  each  gas  species  in  determining  their  cathode  fall  voltage. 
The  density  of  positive  ions  was  always  much  larger  than  the 
density  of  electrons  and  negative  ions  and  thus  determined  the 
slope  -of  the  electric  field  in  accordance  with  Poisson's  equation 
as  well  as  the  voltage  across  the  cathode  fall. 

Conversely,  results  indicated  electron  attachment  in  the 
cathode  fall  region  increased  the  cathode  fall  voltage  and 
length.  Although  the  attachment  coefficient  is  la’rger  than  the 
ionization  coefficient  immediately  in  front  of  the  cathode,  the 
density  of  negative  ions  does  not  contribute  significantly  to  the 
apace  charge  until  the  negative  glow.  The  voltage  drop  and 
length  of  the  cathode  fall  must  increase  in  order  to  maintain  the 
same  net  ionization  as  occurred  in  a  pure  electropositive  gas. 
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